Pacific bluefin tuna, Thunnus orientalis, exhibits a flexible feeding ecology in the Southern California Bight

Pacific bluefin tuna, Thunnus orientalis, migrates from spawning grounds in the western Pacific Ocean to foraging grounds in the California Current System (CCS), where they are thought to specialize on high energy, surface schooling prey. However, there has been substantial variability in estimates of forage availability in the CCS over the past two decades. To examine the foraging ecology of juvenile T. orientalis in the face this variability, we quantified the diet and prey energetics of 963 individuals collected in the Southern California Bight (SCB) from 2008 to 2016. Using classification and regression tree analysis, we observed three sampling periods characterized by distinct prey. In 2008, T. orientalis diet was dominated by midwater lanternfishes and enoploteuthid squids. During 2009–2014, T. orientalis consumed diverse fishes, cephalopods, and crustaceans. Only in 2015–2016 did T. orientalis specialize on relatively high energy, surface schooling prey (e.g. anchovy, pelagic red crab). Despite containing the smallest prey, stomachs collected in 2009–2014 had the highest number of prey and similar total energetic contents to stomachs collected in 2015–2016. We demonstrate that T. orientalis is an opportunistic predator that can exhibit distinct foraging behaviors to exploit diverse forage. Expanding our understanding of T. orientalis foraging ecology will improve our ability to predict its responses to changes in resource availability as well as potential impacts on the fisheries it supports.


Introduction
Pacific bluefin tuna, Thunnus orientalis, migrates from spawning grounds in the western Pacific at ages 0-3 years to forage in the California Current System (CCS) before returning west at ages 3-7+ years [1,2]. These juveniles, typically ranging from~40-175 cm fork length (FL), support commercial and recreational fisheries in both the U.S. and Mexican exclusive economic zones (EEZs), with peak U.S. landings from the Southern California Bight (SCB) [3]. Within the CCS, T. orientalis migrates seasonally from as far south as Baja California, Mexico, in winter to as far north as Washington state in fall [1,4]. Seasonal

Stomach collection, processing, and specimen morphometrics
The research and collection of Thunnus orientalis was conducted under the State of California Department of Fish and Wildlife permit No. SC-12372 issued to the National Oceanic and Atmospheric Administration's Southwest Fisheries Science Center (NOAA SWFSC).
T. orientalis stomachs were collected from 2008 to 2016 between January and October (76% between June and September). Individuals were captured using rod and reel with hook and line by recreational anglers on commercial passenger fishing vessels and private boats. Specimens from at least 121 unique trips were donated by recreational fishers to the Southwest Fisheries Science Center and frozen at -20˚C (S1 Table). Fishing trips operated in the SCB within 150 km of the port of San Diego (Fig 1A). Fork length (FL) and operculum length (OL) in centimeters, as well as capture date and location were recorded for each specimen when available (Fig 1). All fishes were captured during daylight hours, but time of capture was not reported. Capture locations were reported as coordinates (decimal degrees) or the colloquial names of common fishing grounds, which were then converted to approximate coordinates. For specimens where only the head and stomach were donated, FL was estimated from OL [18].
Carcasses were thawed overnight before we removed the stomachs and rinsed contents through a 0.5 mm brass mesh sieve. Stomach contents were stored in 70% ethanol until they could be identified to the lowest possible taxonomic level using published keys [19][20][21][22] and enumerated. No prey items were individually weighed, as most of the prey were at least partially digested and often represented only by hard parts (such as fish bones and cephalopod beaks). Only a single family of hyperiid amphipods was reliably identified (Phronimidae) and was increasingly identified as processor skill increased, so all hyperiid amphipods were grouped by suborder in our analyses (Hyperiidea). Other prey not identified at least to family were combined and considered per broad taxonomic grouping (fishes, cephalopods, and crustaceans). Bait used by fishers (fresh S. sagax and E. mordax >110 mm SL) and stomachs without prey were excluded from all analyses.

Diet analysis
Prey counts were grouped by family and transformed into proportional abundances per stomach prior to diet similarity analyses. This approach reduces the influence of prey with high abundance but low frequency of occurrence and minimizes the effects of predator size on overall prey importance. Proportional prey abundance (p) is calculated by dividing the number of prey per group in an individual stomach by the total number of prey across all prey groups in that same stomach. Only prey families that contributed more than 1% mean proportional abundance (� p) across all stomachs were included in the similarity analyses. Families contributing less than 1% � p were lumped together into an "Other" prey group per broad taxonomic grouping (e.g. "fishes" and "cephalopods").
Tunas have been shown to exhibit diet variability with respect to size, season, capture location, and environmental conditions [12,23,24]. To examine the effects of collection year, month, location, and individual FL on T. orientalis diet, we performed classification and regression tree (CART) analysis in the diet package in R [25]. Collinearity among explanatory variables was examined using the corrplot package [26] (S1 Fig). CART analysis, described in Breiman et al. (1984), is a non-parametric modeling approach with applications to diet analyses extended by Kuhnert et al. (2012). By partitioning data through successive splitting that seeks to minimize an error criterion, the algorithm identifies optimal splitting values for each explanatory variable without requiring a piori knowledge. This allows for empirical partitioning of samples into groups with similar diets (e.g. tuna size classes) as well as non-linear relationships between explanatory variables and diet composition (e.g. variability in tuna size class definitions under different conditions). Additionally, CART analysis can accommodate missing covariate values by splitting on "surrogate" variables that have high correlation to the primary split [27]. These surrogates are also used to calculate the importance of each covariate relative to the selected variable at each split, which allows for recognition of covariates that are good predictors but masked from the final tree. Model performance was assessed using tenfold cross validation and the final model was selected using the "1-SE rule" [28].
Terminal nodes of the CART analysis were used to inform further analyses on groups of T. orientalis with distinct diets. Diet composition was also summarized across terminal nodes of the CART analysis as the percent frequency of occurrence [%FO = (number of stomachs containing prey i /total number of stomachs) � 100] and the percent numerical abundance [%N = (number of prey i /total number of prey) � 100] of each prey type to facilitate comparison with previous studies.
Diet diversity and sample coverage were quantified with diversity accumulation curves using the iNEXT package in R [29]. Diet diversity is represented as the Shannon Index (qD, Hill number of order q = 1) [30]. Diet similarity within and among CART-defined groups was quantified using pairwise permutational multivariate analysis of variance (PERMA-NOVA) in Primer v7 [31]. Homogeneity of multivariate dispersion among groups was quantified using PERMDISP in Primer v7, a multivariate analog to Levene's test for homogeneity of variance [31]. Principal coordinate analysis (PCoA) was performed in the labdsv package in R [32] to visualize the relative contributions of variable means and variances to diet differences among CART groups. Unless otherwise indicated, descriptive statistics are reported as means (± SD).

Prey size and energetics
Lengths were measured for fresh and partially digested prey when possible and reported in mm as standard length (SL) for fishes, mantle length (ML) for cephalopods, and total length (TL) for crustaceans. Whole lengths from partial remains (teleost vertebrae, cephalopod beaks, and crustacean carapaces) were estimated following Glaser et al. [33] using published regressions (S2 Table). All cephalopod MLs were estimated from beaks.
Estimates of total stomach content mass and energetic value required length measurements from all prey. Mean lengths of measured prey were applied hierarchically to unmeasured prey of the same taxon from individuals in: 1) the same stomach; 2) the same sampling year; or 3) the same CART group. Family-level mean lengths were similarly applied to family-level prey IDs. Mantle lengths of unidentified squids for which rostral lengths were directly measured were estimated as the average ML from all beak to body length regressions used for identified cephalopods (S2 Table). For unmeasured, unidentified prey in each group, annual or CART group mean lengths were applied.
To estimate the energetic value of individual prey items, prey lengths were first converted to estimated body mass (g) using published regressions (S2 Table). For cephalopod prey, body mass was estimated directly from beak lengths. The mass of unidentified squids for which rostral lengths were directly measured was estimated as the average mass from all beak to body mass regressions used for identified cephalopods. Body mass was converted to energetic content using published energy density values for prey species following Glaser et al. [34] (S2 Table). Genus-or family-level estimates were used when prey taxonomic resolution was low or species-specific energy estimates were unavailable. For prey not identified at least to family (17% of all prey), energetic values were estimated as annual or CART group-level averages within the relevant taxonomic group and the energetic value of the total reconstituted mass within each stomach was quantified.
Comparisons of total content mass and energetic value among stomachs may be biased if contents represent different amounts of time since the last feeding event. Several studies have shown that bluefin tuna stomach fullness and prey digestion vary with time of day [12,23,35]. However, fullness is a metric based mostly on fresh prey and estimated peak feeding times are highly variable both within and across studies. Although we don't know the exact capture times for any of our specimens, potential bias due to feeding time is reduced by considering the whole mass of both fresh and partially digested prey. Inclusion of reconstituted prey mass limits our temporal bias to the total gut clearance rate, which is also prey dependent and has been estimated at a minimum of 14-20 hours [36,37]. However, the percent of empty stomachs was consistent across sampling years in this study (25 ± 5%, Fig 1B), suggesting a relatively consistent relationship between sampling times and gut clearance.
Our methods could also have overestimated total stomach content mass if mean prey lengths were not representative of unmeasured prey, or we included prey consumed on multiple days. The maximum reported daily rations for bluefin tunas in the size range sampled in this study are~2-3% body mass (BM) [11,12]. To minimize the chances that our methods for assigning masses to unmeasured prey overestimated prey mass or included prey from multiple days, we excluded stomachs with estimated contents > 3% BM from prey mass and energetic analyses.

Comparing prey characteristics across bluefin sizes and CART groups
The size structure of T. orientalis in the SCB was variable across our study period and was reflected in the size structure of the tuna we sampled [3,18]. To examine whether differences in stomach contents across years was impacted by predator size, we quantified the effects of T. orientalis mass on the total mass of prey per stomach as well as individual prey length using generalized additive models (GAMs) in the mgcv package [38]. GAMs were fit using restricted maximum likelihood parameter estimation and smoothing parameters for both models were estimated with Gamma distributions and "log" link functions.
Variability in T. orientalis foraging ecology across our time series was examined by comparing mean prey length, number, and energetic value per stomach among CART groups. None of these comparisons met both basic assumptions of normality and homoscedasticity for parametric mean comparisons, and non-parametric alternatives were used. When data variance was homoscedastic among groups but not normally distributed, a Kruskal-Wallis (KW) rank sum test was performed in R. If the global KW test statistic was significant, differences between group levels were quantified using pairwise Wilcoxon rank sum tests with a Bonferroni adjustment for multiple comparisons [39]. When data were both non-normal and heteroscedastic, a Welch's ANOVA [40] was performed on ranks [41], and pairwise differences between ranks among groups were examined using Games-Howell post hoc tests in the rstatix package [42,43].
The degree of specialization of T. orientalis diet was estimated by examining the prey-specific proportion (analogous to 'prey-specific abundance', [44]) of the most frequently consumed prey in each CART group. For each prey group found in at least 20% of stomachs (prey i ), we quantified the mean proportion and mean number of prey i in stomachs that contained prey i . A prey-specific proportion threshold of 0.5 was used to distinguish specialist (> 0.5) from generalist (< 0.5) foraging behaviors.

Diet summary
Of the 963 T. orientalis stomachs sampled, 721 (75%) contained prey ( Fig 1B). Specimens sampled in 2008-2010 were relatively small (61-93 cm FL, 5.43-10.14 kg), and larger individuals across a broader size range (42-174 cm FL, 1.83-102.54 kg) were sampled during 2011-16 ( Fig  1C). A total of 21,189 prey items were identified across 48 groups (Table 1). Three crustacean groups, seven fish families, and seven cephalopod families contributed >1% � p to T. orientalis diet across all years (Table 1). Of these 16 families, eight were represented by a single species, and four by a single genus. An additional 13 fish families and four cephalopod families contributed <1% � p across all samples and were lumped into "Other fishes" and "Other cephalopods" groups, respectively.

Classification and regression tree analysis
Stomach contents of T. orientalis were partitioned into three groups based on collection year: 2008 (n = 88), 2009-2014 (n = 439), and 2015-2016 (n = 194) (Fig 2A). The pruned classification tree had a cross-validated error rate of 0.225 (SE = 0.02, R 2 = 0.78), and collection year was the most important predictor of diet (relative importance = 1). Although collection month was significantly correlated with year (Pearson's R = -0.61, Figs 1C and S1), it was a poor surrogate for collection year in the CART analysis. Month of collection was largely uninformative (relative importance = 0.01), and neither latitude, longitude, nor specimen fork length were reliable predictors of diet in our dataset (relative importance = 0).  Table).

Diet variability among sampling periods
All CART groups reached 99% sample coverage; only 1% of the total expected family-level prey diversity remained unsampled in our dataset. Although diet diversity was highest in the Table 1. Diet composition of Thunnus orientalis. The percent number (%N), number, percent frequency of occurrence (%FO), and frequency are given for each prey item per CART-defined sampling period. Asterisks indicate prey families that were lumped into an "Other" prey group for diet similarity analyses. Prey ID  2008  2009-14  2015-16  2008  2009-14  2015-16 Cephalopods

Fishes
Carangidae   Table). Variability in diet composition within CART groups was significantly smaller than variability among groups (PERMANOVA, t = 8.87-10.27, p = 0.001). Diet similarity was highest within the 2008 group and lowest within the 2009-2014 group (mean pairwise similarity = 43.06% and 20.60%, respectively). Among group similarity was low for all groups (mean pairwise similarity = 3.13-4.70%). There was also significant variability in dispersion among groups (PERMDISP, Global F = 35.92, p = 0.001, pairwise t = 3.99-9.44, p < 0.005), and PCoA demonstrates that both the mean and variance of diet composition contributed to the observed differences among CART groups (Fig 3B).

Effects of tuna size on stomach contents
Of the 721 T. orientalis with non-empty stomachs, FL and prey % BM could be estimated for 644 specimens. The estimated prey mass to predator mass ratio of 52 stomachs exceeded 3% BM and were excluded from prey mass and energetics analyses (7 stomachs from 2008, 30 from 2009-14, and 15 from 2015-16). The mean mass of prey in the remaining stomachs (n = 592) was 58.3 ± 68.7 g, representing 0.5 ± 0.6% BM, both within reported ranges of bluefin stomach content mass and daily ration for the sizes of tuna sampled [11,12,36]. T. orientalis mass varied among CART groups (Welch's ANOVA on ranks, F = 183.04, p < 0.001, Fig 3C), but was not a useful predictor of estimated stomach content mass among sampling years (GAM, Adj. R 2 = 0.04, Table 2  After accounting for prey group and sampling year, T. orientalis mass was also not a useful predictor of prey length (GAM, Adj. R 2 = 0.22, Table 2   were performed under the assumption that variability in the mass of T. orientalis sampled during our study period was not a significant driver of variability in the quantity, size-structure, or composition of prey in a single stomach. We focused on examining differences in prey characteristics among CART groups, which were represented by distinct taxa.

Prey number and energetics
The total number of prey per stomach was different among CART groups (Welch's ANOVA on ranks, F = 31.60, n = 592, p < 0.001, Fig 5A) Fig 5C). The total energetic values of stomach contents estimated with our methods were within the range of reported daily caloric intake of captive T. orientalis [45].

Feeding behaviors
Using prey-specific proportion to estimate feeding behaviors of T. orientalis, we observed variability in the degree and subject of forage specialization among CART groups (Fig 6). In 2008, T. orientalis specialized on Myctophidae (prey-specific proportion > 0.5), with an average of 15.9 individuals per stomach (prey-specific number). T. orientalis sampled in 2009-2014 exhibited more generalist feeding (prey-specific proportion < 0.5) across a relatively diverse prey assemblage. In 2015-2016 T. orientalis specialized on E. mordax or P. planipes, though the smaller P. planipes were more frequently consumed and in higher numbers (mean preyspecific abundances of 12.7 and 24.5, respectively).

Diet summary
T. orientalis forage on a broad diversity of fishes, squids, and crustaceans in the SCB, including high energy epipelagic prey such as anchovies and sardines, as well as large numbers of mesopelagic prey such as myctophids, hyperiid amphipods, and enoploteuthid squids. While our findings differ from the two previous studies of the T. orientalis diet in the SCB [10,19], enoploteuthid squids are also important prey for smaller T. orientalis (< 30 cm FL) in the western Pacific Ocean [23]. Studies of similar sized bluefin from the Mediterranean (T. thynnus) and southwestern Pacific Ocean (T. maccoyii) reported similarly diverse diets, including mesopelagic species associated with deep-scattering layers [12,35]. Prey items observed in this study were generally smaller than 10 cm, consistent with prey lengths of T. maccoyii in the southwestern Pacific [12] and other tuna species at similar sizes [46]. An increase in predator size often translates into foraging on larger items [47], but we did not observe variability in prey length with increasing T. orientalis mass (Table 2 and S3 Fig). Consumption of relatively small prey has been reported across a broad size range for bluefin and other tunas from multiple ocean basins [11,12,46]. Within the sampling periods defined in this study, T. orientalis ranging in mass from 2 to 103 kg (42-172 cm FL) foraged on similar sized prey (< 10 cm, < 10 g) at similar trophic levels (e.g. [48]). This is consistent with the findings of Madigan et al. (2017) who showed no trend in δ 15 N stable isotope values (a proxy for trophic position) in T. orientalis from~70-200 cm FL in the CCS. Smaller, lower trophic level prey tend to have higher cumulative biomass [49] and may respond more quickly to environmental variability than higher trophic level prey. Targeting diverse, relatively small prey may increase the resilience of T. orientalis in a dynamic system, where forage composition and abundance can change drastically over relatively short timescales.

Thunnus orientalis exhibits a flexible foraging ecology in the Southern California Bight
T. orientalis in this study were classified into three sampling periods based on prey taxonomic composition, with limited overlap in the dominant prey species. Additionally, prey consumed during these three sampling periods also had different size structures and behaviors. In 2008, T. orientalis fed mostly on medium sized (3-5 cm) myctophids and enoploteuthid squids, abundant constituents of the daytime deep scattering layer in the southern CCS (~300-500 m depth) [50][51][52] that migrate into near-surface waters at night [53,54]. Archival tag data and diet analyses (where time of feeding was inferred from prey digestion state) suggest that bluefin tunas can feed during day and night [12,35]. Although we do not know at what time of day the prey in this study were consumed, both day-and nighttime depth distributions of myctophids and enoploteuthid squids are within the diving range of T. orientalis (~0 to > 450 m) [5,10].
In 2015-2016, T. orientalis fed heavily on P. planipes and E. mordax (3-8 cm), both of which can form dense schools in near-surface waters [55,56]. These prey were often found in large numbers in stomachs containing little else (Fig 6), suggesting T. orientalis exploits such aggregations. Specialization on abundant, near-surface schooling prey most closely resembles the findings of Pinkas (1971, E. mordax) and Madigan et al. (2015, S. sagax). However, this strategy is not limited to epipelagic prey as T. orientalis also appeared to specialize on myctophids in 2008 when using a prey-specific proportion threshold of 0.5 (Fig 6, [10,44]).
Prey consumed in 2009-2014 were the smallest (1-5 cm), most taxonomically diverse, and were consumed in the highest numbers. It is unlikely the increased diversity was simply due to a longer sampling period as each individual stomach had relatively high diversity (Fig 3A and  S3 Table). The most abundant forage species also had diverse depth habitats, including the epipelagic Pacific jack mackerel, (Trachurus symnmetricus) [57], and the more broadly mesopelagic gonatid squids (Gonatidae) [53] and hyperiid amphipods (Hyperiidea) [58]. In contrast to the periods of prey specialization observed in this and previous studies, T. orientalis could be described as a generalist predator during 2009-2014, feeding throughout a broader swath of the water column in the SCB. This generalist foraging has also been observed in T. thynnus and T. maccoyii [12,35]. Young et al (1997) found that T. maccoyii in nearshore and offshore habitats fed on diverse, but distinct diets that were correlated to differences in prey availability between habitats. While we were unable to directly compare diets to estimates of forage availability, each sampling period was characterized by prey with different depth habitats and schooling behaviors. Thus, it is likely that T. orientalis employs different feeding behaviors to opportunistically exploit these ecologically distinct forage.

Diet variability coincides with responses of forage to climatic variability in the CCS
Although we were not able to quantify the effects of environmental conditions on diet variability in our dataset, temporal variability in T. orientalis foraging ecology coincides with major oceanographic and biological variability in the SCB. During the late 2007 to early 2009 La Niña event, the SCB was characterized by relatively cool sea surface temperatures, accompanied by low to neutral anomalies of phytoplankton and zooplankton biomass [59]. At the same time, relatively high abundances of larval myctophids were also observed [17], suggesting a large spawning biomass that could have been exploited by T. orientalis. The same surveys also noted low abundances of larval anchovy and sardine in 2008. ENSO neutral conditions prevailed from 2011-2014 following brief El Niño and La Niña events in 2009-2010 and 2010-2011, respectively. Consistently positive chlorophyll-a anomalies and variable water column temperatures during this period supported an order of magnitude variability in zooplankton biomass [60] and diverse larval fish assemblages, from which anchovy and sardine were notably absent [61]. Increased sea surface temperatures and stratification in the SCB during a 2014-2015 marine heat wave and an El Niño event in 2015-2016 were associated with relatively high abundances of anchovy [61]. During these warming events, the diets of multiple predators became dominated by anchovy (e.g. California sea lions (Zalophus californianus), rhinoceros auklet (Cerorhinca monocerata), [16]) and pelagic red crabs (yellowtail jack (Seriola dorsalis, pers. communication and observation), reflecting concerted food web changes across the CCS.
While these shifts were apparent in the diets of bluefin, little is known about the relationship between T. orientalis diet and metrics of forage availability in the CCS, the most extensive of which are derived from trawl surveys (e.g. CalCOFI and NOAA's Coastal Pelagic Survey [16,17]). Direct comparisons of predator diets to these trawl surveys are an important step in determining how well trawl estimates of forage availability reflect forage use by predators. For managed forage species (e.g. anchovy, sardine, rockfishes (Sebastes spp.)), predator diet monitoring could refine estimates of natural mortality under different environmental conditions. Ultimately, it will be important to continue diet analyses across seasons and climatic regimes (e.g. El Niño, La Niña) to more directly link predator diets to forage availability and oceanography variability.

Feeding on diverse forage provides similar meal energetic value
Among sampling periods, there was limited variability in estimated meal energetic value (kJ per stomach, Fig 5C) despite differences in prey composition, sizes, and numbers (Figs 2, 4 and 5A). There was no difference in the total energetic value or number of prey per stomach when T. orientalis was specializing on myctophids in 2008 or P. planipes and E. mordax in 2015-2016 (Fig 6). However, the mean prey-specific number of E. mordax (12.7) was lower than Myctophidae (15.9), and the standard deviations of prey-specific numbers were higher in 2015-2016 than 2008. Myctophid aggregations have lower average densities and are larger on average (0.2-200 individuals m -3 , 100s of m diameter by 10s of m thick [62][63][64]) than schools of E. mordax (0.5-500 individuals m -3 ,~10s of m diameter x 20 m thick) [56,65]. Myctophids may represent more consistent patches of prey across larger areas than E. mordax or P. planipes, which could explain their higher frequency of occurrence and lower variance in prey specific number (Fig 6) [56,66,67]. Stomachs collected in 2009-2014 contained prey with relatively low individual energetic values but had similar total energetic values to stomachs collected in 2015-2016 due to the high numbers of prey (Fig 6). Daily caloric intake is determined by the identity and number of prey consumed, but net energetic gains are tied to energy expenditure during search and capture [68].
T. orientalis exhibits diverse diving behaviors across seasons and habitats that have been tied to foraging [6, 7] and likely require different amounts of energy. Foraging near the surface on relatively high-density prey may have lower energetic costs than foraging throughout the water column on more diffusely distributed prey [69]. Additionally, the metabolic rate of juvenile T. orientalis may be elevated while foraging in relatively cool water at depth, outside of their metabolic thermal minimum zone of 15-20˚C [70]. Thus, prey consumed in 2009-2014 may have been more costly to capture than red crabs and anchovies in 2015-2016 and resulted in lower net energetic gains or growth rates. Madigan et al. [71] reported reduced growth rates in T. orientalis that fed higher proportions of midwater prey, although the influence of tuna size on this observation remains unclear. Future diet and archival tagging studies that can more directly link feeding and expenditures to predator condition or growth rate will elucidate the net value of foraging on these distinct prey groups.

Implications for our understanding of bluefin tuna ecology and fisheries
Globally, bluefin tunas are often associated with energetically dense, schooling forage [11,12,19]. The distribution of forage fishes impacts the migration and density of tuna schools, and the availability of high energy forage is positively correlated with tuna condition and landings [6,8,13,72]. Although it is not clear how prey switching affected its horizontal distribution or abundance in the SCB, daytime depth distributions of T. orientalis significantly impact its availability to recreational fishers, which largely operate during daytime. When considering the daytime depths of dominant prey across the sampling periods defined in this study, the daytime depth of T. orientalis likely shoaled throughout our study period (Fig 7). Recreational fishers had a relatively difficult time hooking T. orientalis in 2008-2014, when fish were only found at depths from 50 to 150 meters, compared to 2015-2016 when fish were consistently observed in surface waters jumping and boiling on bait balls (pers. communication and observation). Understanding the foraging ecology of T. orientalis can help to explain differential effort required to capture tunas.

Conclusions
We present the most extensive time series of T. orientalis feeding ecology to date and describe a more flexible diet than previously observed in the CCS. T. orientalis is an opportunistic forager that exhibits periods of specialization and more generalist feeding, demonstrating the importance of multi-year diet studies to capture the dynamics of foraging ecology. This flexible feeding ecology has likely contributed to the persistence of T. orientalis in the CCS in the face of significant ecological variability over the past two decades. Continued diet monitoring that explicitly links diet to oceanographic conditions will improve our understanding of the drivers of variability in T. orientalis foraging ecology and may increase our capacity to predict changes in its distribution and availability to fishers.  Table 2. Each panel shows the relationship between a covariate and the contribution of the smoother for that covariate to the model's fitted values ("s(x)"). Grey shading indicates 95% confidence intervals about the estimate for each covariate, and rug plots indicate covariate observations. (TIF)  Table 2. Panels (a) and (b) show the relationship between a covariate and the contribution of the smoother for that covariate to the model's fitted values ("s(x)"). Grey shading indicates 95% confidence intervals about the estimate for each covariate, and rug plots indicate covariate observations. In panel (c), the effect of prey group on prey length ("f (x)") is given in reference to the mean length of cephalopods and error bars indicate 95% confidence intervals about each group mean. (TIF) S1 Table. Table. Mean proportional abundance (� p) and percent frequency of occurrence (%FO) of prey groups included in the CART analysis for each terminal node ("CART group"). "Other" groups include prey identified at least to family that represented less than 1% (� p) across the dataset (indicated by asterisks in Table 1). (DOCX) S4 Table. Sample coverage and diversity estimates for Thunnus orientalis stomachs in each CART group. N = the total number of stomachs per year group and t = the sample size at 99% sample coverage. Diversity estimates using the Shannon Index are given at 99% sample coverage (qD ± 95% confidence interval).